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00 Abstract. We present a new determination of the mass content of the Sculptor dwarf spheroidal 

^^ galaxy, based on a novel approach which takes into account the two distinct stellar populations 

^^ . present in this galaxy. This method helps to partially break the well-known mass-anisotropy 

vN ■ degeneracy present in the modelling of pressure-supported stellar systems. 
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1. Introduction 



C^ The determination of the mass content of galaxies is a fundamental step for our un- 

Q derstanding of their formation and evolution. This is particularly true for small systems 

5-H ' such as dwarf spheroidal galaxies (dSphs) , whose total mass can determine their destiny 
t^ , as being fragile objects, strongly perturbed by supernovae explosions and/or interactions 
I ^i ' with the environment, or mostly undisturbed objects. 

Part of the interest in determining the masses of dSphs has surely been driven by the 
^~~l ■ fact that these galaxies, with mass-to-light (M/L) ratios up to 100s, appear to be the most 
^ , dark matter (DM) dominated objects known to-date. This makes them potentially good 
^^ testing grounds for different DM theories of galaxy formation; however at the moment 

^fs it is still unclear if dSphs inhabit cusped DM profiles, as predicted by cold DM theories 

fSj ' of galaxy formation, or cored profiles, typical of warmer kinds of DM. This distinction is 

r I , often hampered by the degeneracies which are intrinsic to the method generally used to 

f~^ ■ determine the mass of pressure-supported stellar systems, i.e. a Jeans analysis of the line- 
OO , of-sight (l.o.s.) velocity dispersion obtained considering all stars as a single component 
^D ■ embedded in an extended DM halo. This analysis is subject to the well known degeneracy 
^ , between the mass distribution and the orbital motions of the individual stars in the 
system (mass-anisotropy degeneracy). Recent observations have shown that some dSphs 
host multiple stellar populations which have distinct spatial distribution, mctallicity and 
j^ ' kinematics. The presence of multiple stellar populations requires a modification in the 
way these systems are dynamically modelled to derive their mass content. 

Here we describe the results of our determination of the mass content of the Sculptor 
(Scl) dSph by adopting a more detailed kinematic modelling than the traditional one, 
i.e. by considering Scl as two- (stellar) components embedded in an extended DM halo. 
We first briefly describe and discuss the observational properties of Scl derived from our 
work (Tolstoy et al. 2004; Battagha 2007; Battagha et al. 2008). We then move to the 
results of the two-component modelling; we show a clear example of the mass-anisotropy 
present when modelling Scl as a one-stellar component system, and finally discuss our 
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results. For a more detailed description we refer to Battaglia (2007) and Battaglia et al. 
(2008). 



2. Observed properties of the Sculptor dSph 

We acquired extended photometry from the ESO/2.2m WFI and VLT/FLAMES spec- 
tra (R~6500) in the Call triplet region for hundreds of individual red giant branch 
(RGB) stars in Scl, out to its nominal tidal radius (Tolstoy et al. 2004: 308 probable 
members; Battaglia et al. 2008: 470 probable members). We used these data to study 
the large scale kinematics and metallicity properties of Scl, and the spatial distribution 
of its stellar populations. 

The combined information from the accurate line-of-sight velocities (±2 km s^^ ) and 
metallicities (±0.15 dex) of RGB stars allowed us to unveil the presence of two distinct 
stellar populations with different spatial distribution, metallicity and kinematics: the 
metal rich (MR) stars ([Fe/H]> —1.5) are more centrally concentrated and show colder 
kinematics than the metal poor (MP) stars ([Fe/H]< —1.7), as can be seen from the line- 
of-sight (l.o.s.) velocity dispersion profiles in Fig. [TJ The MR and MP RGB stars appear 
to trace, respectively, the red and blue horizontal branch (RHB, BHB) stars. The surface 
number density profile of RGB stars, representative of the overall stellar population of 
Scl, is well approximated by a two-component fit, where each component is given by the 
rescaled best-fitting profile derived from our photometry for the RHB and BHB stars. In 
contrast, a single profile (be it a Plummer, Sersic or King) does not appear to be a good 
approximation for the surface number density of RGB stars. Multiple stellar populations 
are known to exist also in other dSphs, such as Fornax (Battaglia et al. 2006) and Canes 
Venatici (Ibata et al. 2006), and might be a common feature to this class of objects. 

The large spatial coverage and statistics of this spectroscopic data-set has allowed 
the discovery of another interesting feature, i.e. a velocity gradient of 7.6 km s~^ deg~^ 
in the Galactic standard of rest (GSR) IjJ along the projected major axis of Scl. This 
gradient is likely to due intrinsic rotation, and not to tidal disruption (see Battaglia 2007 
and Battaglia et al. 2008 for discussion on this point), and this is the first time that 
statistically significant rotation is found in a dSph. 

Among the other Milky Way (MW) dSphs, a significant velocity gradient is detected 
in Carina (Munoz et al. 2006) and a mild one in Leo I (Mateo et al. 2008), but these are 
interpreted as the result of tidal interaction with the MW. However, velocity gradients 
are not exclusive of the dSphs satellites of larger galaxies, but are present also in some of 
the isolated Local Group dSphs, such as Cetus (Lewis et al. 2007) and Tucana (Fraternali 
et al., in preparation), which are unlikely to have ever interacted with either the MW or 
M31 and therefore suffered tidal disturbance. This might point to rotation as an intrinsic 
property of this class of objects and provide insights into the mechanisms that shaped 
their evolution. This might for instance favour scenarios in which the progenitors of dSphs 
were rotating systems which were transformed into prevalently pressure supported object 
by tidal stirring from their host galaxy (e.g., Mayer et al. 2001). It would therefore be 
interesting to compare the frequency of occurrence of rotation among dSphs and other 
generally more isolated systems such as dwarf irregular galaxies and transition types. It 
is possible that more velocity gradients will be discovered in MW dSphs now that larger 
and more spatially extended velocity surveys are becoming available. 

In the analysis presented here we always use rotation-subtracted GSR velocities. 

f We use the velocities in the GSR frame to avoid spurious gradients introduced by the 
component of the Sun and Local Standard of Rest motion's along the l.o.s. to Scl. 
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3. The Jeans analysis 

DSphs generally have low ellipticities and are presumably prevalently pressure sup- 
ported, and therefore are commonly considered as spherical and with no streaming mo- 
tionqj. Then, if the system is in dynamical equilibrium, the l.o.s. velocity dispersion 
predicted by the Jeans equation is (from Binney & Mamon 1982) 



LiR)-^rjj^, I '-i^=fi^-P-a)dr (1) 



where R is the projected radius (on the sky) and r is the 3D radius. The l.o.s. velocity 
dispersion depends on: the mass surface density S*(_R) and mass density p*(r) of the 
tracer, which in our case are the MR and the MP RGB stars; the tracer velocity anisotropy 
/?, defined as /3 = 1 — a^ja^, and its radial velocity dispersion ar,*^ which depends on 
the total mass distribution (for dSphs the contribution due to stellar mass is negligible). 
A unique solution for the mass profile can be determined from this equation if both 
the mass distribution of the tracer and ii{f) are known, although this solution is not 
guaranteed to produce a phase-space distribution function that is positive everywhere. 
However, proper motions for individual stars in the system would be necessary to derive 
/3(r) and this is not feasible yet. This causes the well-known mass-anisotropy degeneracy, 
which prevents a distinction amongst different DM profiles. Examples of this are present 
in the literature: the observed l.o.s. dispersion profiles of dSphs are both compatible with 
cored DM models, assuming that /3 = (e.g., Gilmore et al. 2007), and with cusped 
profiles, allowing for mildly tangential, constant with radius, /3 (e.g. Walker et al. 2007). 
It is therefore advisable to explore several hypotheses for the behaviour of /3(r) and for 
the mass distribution when carrying out a Jeans analysis (see below). 

In the following we explore the possibility that DM follows a cored (pseudo-isothermal 
sphere) or a cuspy distribution (NFW profile); we also explore two hypotheses for the 
behaviour of /3(r): constant with radius, and following an Osipkov-Merritt (OM) model 
(isotropic in the central parts and radial in the outskirts). The spatial distributions of 
the tracers are derived from our photometric data. 

Two- components modelling Here we show the results of modelling Scl as a two-(stellar) 
component system embedded in a DM halo. Eq. [T] is applied separately to the MR and 
MP stars, assuming that they move in the same DM potential. We explore a range 
of core radii Vc for the pseudo-isothermal sphere (re = 0.001,0.05,0.1, 0.5, 1 kpc) and of 
concentrations c for the NFW profile (c = 20, 25, 30, 35). By fixing these, each mass model 
has two free parameters left: the anisotropy and the DM halo mass (enclosed within the 
last measured point for the isothermal sphere, at 1.3 deg = 1.8 kpc assuming a distance 
to Scl of 79 kpc, see Mateo 1998; and the virial mass in the case of the NFW model). We 
compute a x^ for the MR and MP components separately (xmr ^^'^ Xmpj respectively) 
by comparing the various models to the data. The best-fit is obtained by minimising the 

sumx^ =Xmr + Xmp- 

Figure [T^,b shows the models with constant anisotropy, which can be excluded as 
in this case neither a cored or a cusped model is a good representation of the data 
(Xmin "^ I''')- The situation is different in the hypothesis of an OM velocity anisotropy 
as shown in Fig. [T};,d). We find that a pseudo-isothermal sphere with a relatively large 
core {re =0.5 kpc, M(< 1.8 kpc) = 3.4 ± 0.7 x 10^ M©) gives an excellent description 
of the data (Xmin = ^.9). Also an NFW model is statistically consistent with the data 

X We checked that the assumptions of sphericity and absence of streaming motions have a 
negligible effect on the results. 
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Figure 1. L.o.s. velocity dispersion profile (squares with errorbars), from rotation-subtracted 
GSR velocities, for the MR (left) and MP (right) RGB stars in Scl. The lines show the best-fitting 
pseudo-isothermal sphere (solid) and NFW model (dashed) in the hypothesis of /3 = const (a,b) 
and /3 — /3om (c,d). The models have reduced x^ ~ 1-6 (a,b), 0.6 (c) and 1 (d). 



(c = 20, virial mass M„ = 2.2;o° x 10^ Mq, xLin = 10.8), but tends to over-predict the 
central values of the MR velocity dispersion. The mass within the last measured point 
(~ 2.41q 7 X 10^ Mq) is consistent with the mass predicted by the best-fitting isothermal 
sphere model. 

One- component modelling Figure [2] shows a clear example of the mass-anisotropy de- 
generacy when modelling Scl with the traditional one-component analysis. Just by chang- 
ing the assumption on /3 (=const. or — /3om), cored profiles with very different core radii 
and M(< 1.8 kpc) give excellent fit to the data. An excellent fit is also an NFW profile, 
in the hypothesis oi P — const. These models are practically indistinguishable and all 
good (x^ ^ 8, reduced x^ ~ l-l)- No DM model can be favoured nor hypotheses on /3. 



4. Discussion and conclusions 

The combined information of velocity and metallicity from a large sample of FLAMES 
spectra of individual stars in Scl has allowed us to separate the different kinematics of 
the two stellar populations in Scl, and to carry out a more detailed kinematic modelling 
which allows us to partially break the mass-anisotropy degeneracy. 

Under the hypotheses explored for the behaviour of the velocity anisotropy f3{r), the 
two-component modelling indicates that a cored profile is slightly favoured over a cusped 
profile, as the latter tends to over- predict the central values (at i? < 0.1 kpc) of the 
MR velocity dispersion. It could be argued that the relatively large central dispersion 
predicted for the MR stars from the cusped profile could be lowered by allowing the 
velocity anisotropy of the MR population to be tangential at i? < 0.1 kpc. However it 
should be noted that in order to reproduce the rapid decline of the MR velocity dispersion 
profile, the MR (3{r) needs to become very radial already on the scale of 0.3-0.4 kpc, given 
the observed spatial distribution of MR stars. Observational determinations of /3, either 
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Figure 2. a) Observed l.o.s. velocity dispersion profile, from rotation-subtracted GSR velocities, 
for all RGB stars in Scl (squares with error-bars), with overlaid the best-fits for a cored DM 
distribution (solid line: /? — /3om and r^ =0.5 kpc; dotted line: /3 = const. and r^ =0.05 kpc), 
and an NFW DM distribution (dashed line: (3 = const, and c = 35). Panel b) and c) show the 
velocity anisotropy and mass distribution corresponding to the models in panel a). Note that 
the model with r^ =0.05 kpc is almost equivalent to a cusp. 



from direct measurement or from distribution functions, would be the most desirable 
(and definitive) solution to this issue. 

Strigari et al. (2007) showed that, independently of /3(r) and the total mass distribu- 
tion of dSphs, the mass enclosed within 0.6 kpc (Mq.s) can be derived more accurately 
than at any other distance from the centre. However, given that the evolution of dSphs 
is most likely to be governed by their total mass and that dynamical analyses provide 
only the mass enclosed within the last measured point, it is important to have reliable 
determinations of the mass enclosed as farther out as possible. Our analysis shows that, 
independently of the adopted mass profile, the best fits from our two-component mod- 
elling yield consistent values for the mass enclosed within the last measured point, at 
l.Skpc. 

Finally we would like to comment on the issue that dwarf galaxies might be embed- 
ded in haloes of similar mass (~ 4 x 10^ M©) as has been claimed in the literature. 
Figure [3fleft) shows a revisited version of the mass enclosed within the last measured 
point versus absolute magnitude plot. The masses plotted are the ones derived in Walker 
et al. (2007) for a sample of dSphs, because they are derived from the most extended 
velocity samples and from a homogeneous analysis. For Scl we use our own determina- 
tion, which is the most accurate so far. The figure shows that the DM mass of Scl is now 
much larger than the typical DM mass derived so far for the other dSphs. Also Fornax 
deviates considerably from the relation. In the plot shown here, the masses of Fornax and 
Sextans were derived from data extending out to ~0.6 and 0.3 tidal radii, respectively. 
Given that Fornax and Sextans have been recently mapped much farther out (Fornax 
out to ~ 1 tidal radius: Battaglia et al. 2006; Sextans out to ^ 0.8 tidal radius: Battaglia 
et al. in preparation), their revised masses are likely to increase quite significantly, and 
introduce even larger dispersion in this plot. In some cases, Mo.e has been used in this 
kind of analysis, instead than the mass enclosed within the last measured point. However, 
Fig. [Sjright) shows that Mo.e is very insensitive to the total mass of the dSph. In fact, 
NFW haloes of virial masses between 2-40 xlO^ M0(i.e. in the range found by Walker 
et al. 2007 to fit their samples of dSphs) give remarkably similar values of Mq.g, probably 
because at such a small distance there is not much sensitivity to the total mass. Therefore 
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Figure 3. Left: Mass within tiie last measured point (rmax) versus absolute V magnitude, M^, 
for MW dSphs. The masses are from Walker et al. 2007 (squares), except for Scl which is from 
our two-components modelling (asterisk); M^ is from Irwin & Hatzidimitriou (1995). The arrows 
indicate that the masses measured for Sextans and Fornax are likely to increase given that these 
objects have been surveyed just to 0.3 and 0.6 tidal radii. The solid line indicates amass of 4x10^ 
M0. Right: Mass enclosed within r for NFW halos of virial masses 2,4,8,20,40x10** M0(solid 
lines), i.e. the range found by Walker et al. (2007) to fit the observed l.o.s. velocity dispersion 
profile of their sample of 7 dSphs. Concentrations are derived using the formula in Koch et al. 
(2007), extrapolated from Jing & Suto (2000). The dashed lines enclose masses between 2-7x 10 
Mq; the dotted line is placed at a radius of 0.6 kpc. 

there is no indication of a common mass scale amongst the dSphs in the Local Group, 
while perhaps one may draw the conclusion of the existence of a minimum mass for the 
formation of dSphs. 



References 

Battaglia, G., et al. 2006, A&A, 459, 423 
Battagha, G. 2007, PhD thesis, 



Univ. 



Groningen, The Netherlands, 



http://irs.ub.rug.nl/ppn/304002712 
Battaglia, G., et al. 2008, ApJ (Letters), 681, L13 
Binney, J. & Mamon, G. A. 1982, MNRAS, 200, 361 
Gilmore, G., et al. 2007, ApJ, 663, 948 

Ibata, R., Chapman, S., Irwin, M., Lewis, G., & Martin, N. 2006, MNRAS (Letters), 373, L70 
Irwin, M. J., & Hatzidimitriou, D. 1995, MNRAS, 277, 1354 
Jing, Y. R, & Suto, Y. 2000, ApJ (Letters), 529, L69 
Koch, A. et al. 2007, ApJ, 657, 241 
Lewis et al. 2007, MNRAS, 1364, 1370 
Mateo, M. L. 1998, ARAA, 36, 435 

Mateo, M., Olszewski, E. W., Walker, M. G. 2008, ApJ, 675, 201 
Mayer, L., et al. 2001, ApJ, 559, 754 
Mufioz, R. R. et al. 2006, ApJ, 649, 201 
Strigari, L. E., et al. 2007, ApJ, 669, 676 
Tolstoy, E., et al. 2004, ApJ (Letters), 617, L119 
Walker, M. G., et al. 2007, ApJ (Letters), 667, L53 



